CLAIMS 

1. Method of determining the,, velocity V and 
anellipticity r| parameters for processing seismic traces 
in a common midpoint. (CMP) * gather including an anelliptic 
NMO correction, comprising: . • J«. - : 
5 . a preliminary step' to i define a- plurality of nodes 

(dtn, t 0 ) , the said nodes- being indicative of parameters 
dtn. and x 0 v representing the NMO correction for - the . 
maximum - offset"': and;:; the "t 'zero 'off set '^travel tittie^ in: 
hyperbolic coordinates f v the".' said preliminary ^ being 
10 followed by . ^ . - ■ ■ 

- for each node (dtn/ t 0 ) defined in the preliminary 
step, the following steps: 

- for static NMO correction of traces in the CMP 
gather as a function of .the values of the said parameters 

15 dtn, and t 0 at the .node considered /: and . .. 

- for calculating the,,, semblance function 
associated 1 with the said NMO correction K for the node 

'*' "* considered ;V " aLfrd^v. r- ; :\ ZZZ>. .. f )\ ' : y-^-y.i ") *,. Z* \ Z- ZZ-Z' [ Z - 

- for each picked time t 0 , a step including 

20 determination of the maximum semblance node (dtn(t 0 ) , 

. : - - '.' : .0 : - . 'r.Z^:ZZ:ZZ^. ■ "Z ^ZZ^ : ^M^Z^^Z-t\ \iZ: .• ' r • ■ 

TO (t 0 ) ) / : - . • • ; " / , . 

- and a final.; step, tb^ ^ 

parameters/, so . as ... tc^ obtain y the. ...Velocity y V(to). and 

anellepticity T) ( 1 0 ) -laws*. \: 'V r Y '':ZZ\''.. 'Z-y^. 

25 2. Method according to claim 1, wherein- the nodes 

are defined during the preliminary step in an analysis 
volume (dtn, x 0l t 0 ) determined by minimum and maximum 
values respectively [dtrw, dtn^xl fromim Tomaxl and [tonin, 
tomax] of the dtn, x 0 and t 0 parameters . 
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3. Method according to claim 2, wherein, during the 
preliminary step, a corridor [dtiw (t 0 ) , dtn^fto)], 
[Tominfto), Tomax (to) ] for changing dtn and To parameters is 
delimited inside the analysis volume as a function of 
plausible velocity V and anellipticity n. values, the 
nodes (dtn, t 0 ) defined for. applying ! 'the""NMO correction 
being then located along the. corridor thus delimited. 

4. Method ' according v to any "one of the preceding 
claims > further comprising, for each node (dtn, x 0 ) , a 
stacking step of the corrected seismic traces, following 
the semblance function calculation step, i . 

5. Method according to claim 4, wherein the stacking 
of corrected traces is done using only near offset 
traces. 

6. Method according to any one of. claims 4 or 5, 
further comprising for each picked time," and following 
the step for determining the maximum semblance node , a 
step of checking that values dtn and t 0 of the maximum 
semblance node correspond to a stacking extreme value for 
the same, values dtn and t 0 '; •'• 

7. Method, according to any one of the preceding 
claims, further comprising a step of selecting and 
adjusting the pickings obtained, following the step 
implemented for determining the maximum semblance node 
(dtn (t 0 ), T 0 (t 0 )) for each picked time t 0 , before the 
conversion step. : "' 

8. Method according to claim 7, wherein the said 
step' of selecting and adjusting the pickings comprises a 
step of only retaining pickings dtn and t 0 for which time 
to the highest semblance pickings is greater than a 
predefined value. 



9. Method according to claim 8, wherein the said 
step of selecting and adjusting the pickings also 
comprises a step for adjusting the retained pickings dtn 
and t 0 by parabolic interpolations using values about the 

5 said picked values. 

10. Method according to claim ?, wherein, the said 
step of selecting and adjustingt pickings also comprises 
a step of eliminating retained and adjusted pickings dtn 
and t 0 when it is impossible to calculate the Dix interval 

10 velocities, between the picking considered^ and, higher, 
semblance pickings. 

11. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic traces is. an 
NMO correction process implementing a static correction 

15 CORRnmo- 

12. Method according to claim 11, wherein, during 
the preliminary . step, the NMO corrections, CORRnmo are 
calculated, for all, nodes (dtn, t 0 ) included in the 

• ^aly^ volume and all offsets 'of - processed ^seismic 
20 ' traces . . '. ; : ; ; " ; \/v ' ^- J -\.;\'- K' k: '^' ." ■ Vv'-/ 

13 . Method according to claim 12, wherein the NMO 
correction carried out for each node (dtn, t 0 ) consists of 
app lying NMO . correct ions CORRimo calculated" .during t he 
preliminary step . ,'■ \ v : ,- , . .. Lt?\. : 

•25 l 4 ' : Method according to any one of claims. 11. to 13, 

wherein for a given (dtn, t 0 ) pair, the static NMO 
correction CORR» of .. a seismic trace with , offset x is 
carried out according to the following equation: 

CORWx)- - " To. + U dtn(dtn-f2T 0 ) x2 ^ ^ 

V X max 

3 0 represents the maximum Offset in the CMP gather. 
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15. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic traces is ' a 
PSTM migration using a static NMO correction CORR PSTM . 

16. Method according : to claim 15, wherein, ■ during 
5 the preliminary step, ' the NMO corrections' CORR PSTM are 

'calculated for all nodes, (dtn., and .' x 0 ) included, in the 
analysis volume and all migration off sets inside the 
migration aperture. 1 ^' - v ■ . „ ? 

v 17. Method according to claim 16, wherein tlie NMO 
10 correction step carried out for; e ; ach node (dtn : and x 0 ) , 
comprises, for each offset class, application. of the said 
NMO corrections CORR PSTM , calculated . during the. 
preliminary step on all midpoints inside the migration 
aperture. 

15 18. Method according to claim 17, wherein the NMO 

correction step carried out for: each node (dtn and t 0 ) 
comprises, for each offset; class, ': the- stack : of the 
corrected midpoints following application of the. said NMO. 
corrections CORR PSTM . „ - ' •"-> - „ ' \ * v « 

20 . . 19 . V ; Me t ho d a c c or d i ng to /ah^ • 
wherein, for a given pair (dtn and ? t 0 J > the static NMO 
correction CORR PSTM is carried out : according to the 
following equation : \ f ^-X>,! 

Cor r psTM (x) = -To + fiF, dtn(din + tejs~x + W : r '+ - ^ dhi(dtn%2r 0 j(x-x-h^. 

25 ' ' > 'whete: ■ -v^'-WV. 

- . represents the., coordinates of the midpoints, • 

- x - Xn represents the migration aperture PSTM, 

- h is the half source - receiver offset, 

- xmax is the maximum offset and aperture of the 
3 0 migration. 
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20. Method according to any one of claims 14 or 19, 
wherein, during the final conversion step, the parameters 
dtn .(to) and (xO) are converted' to the velocity law V(t0) 
according to the following equation: 

x- 



V = 



max 



jdtn(dtn + 2T 0 ) — 



V 



21. Method according to one of claims: 14 • or 19, 
wherein, during the final conversion , step, the parameter 
TO (to) ^ is converted to, ; the ane 1 leptic ity . r\ (tp, ) : 1 aw 



according to t| = i 



10 22. Method according' to claims 20 and 21, wherein 

parameter dtn is defined with respect to the velocity V 
and anellepticity r| . according to the following equation: 




l Q , A max 



1 + 8^):/ (l + 8r|)V 2 : 

15 v 23. Method according to claim 21, ; wherein parameter 
, xO is defined according/ to; anellepticity > : . t] according, to 
the following equation: • - : • vr - v • 



To = 



1 + 8t| 



. 24. Method . of. characterising a velocity field for, 

20 processing seismic data, using 

at common midpoint , wherein >• for each travel time to. for 
a zero offset, a set of parameters dtn and t0 is defined, 
representing the NMO correction for maximum offset, and 
the zero offset travel time respectively, in hyperbolic 

25 coordinates. 
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